hep-ph/9405366 



ttN Scattering in a Large- N c cr-Model 

Michael P. MattisQ and Richard R. Silbar 

Theoretical Division, Los Alamos National Laboratory 
Los Alamos, NM 87545 USA 
(May 1994) 



Abstract 



We review the large- N c approach to meson-baryon scattering, including 
recent interesting developments. We then study ttN scattering in a particular 
variant of the linear cr-model, in which the couplings of the a and tt mesons 
to the nucleon are echoed by couplings to the entire tower of I = J baryons 
(including the A) as dictated by large-iV c group theory. We sum the com- 
plete set of multi-loop meson-exchange irN — * ttN and irN — * aN Feynman 
diagrams, to leading order in 1/N C . The key idea, reviewed in detail, is that 
large-iV c allows the approximation of loop graphs by tree graphs, so long as 
the loops contain at least one baryon leg; trees, in turn, can be summed by 
solving classical equations of motion. We exhibit the resulting partial- wave S- 
matrix and the rich nucleon and A resonance spectrum of this simple model, 
comparing not only to experiment but also to irN scattering in the Skyrme 
model. The moral is that much of the detailed structure of the meson-baryon 
S"-matrix which hitherto has been uncovered only with skyrmion methods, 
can also be described by models with explicit baryon fields, thanks to the 
1/N C expansion. 
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I. OVERVIEW OF THE LARGE- iVc APPROACH TO MESON-BARYON 

SCATTERING 



It is well known (l|-|6j that QCD simplifies greatly in the limit N c — > oo, N c being the 
number of colors. Not surprisingly, the large- N c limit has likewise proved to be very useful 
in studying effective low-energy hadron Lagrangians for the Strong Interactions. Broadly 
speaking, such effective theories fall into two categories. On the one hand, there is the 
straightforward Feynman diagrammatic approach in which mesons and baryons are each 
treated as explicit dynamical fields, while on the other hand, there is the more economical 
skyrmion picture J/]|| in which baryons are viewed as solitons constructed from the meson 
degrees of freedom. Since both these approaches purport to describe the low-energy Strong 
Interactions, it follows that if they are sensible, they should be equivalent to one another. 
Furthermore, this equivalence must hold order by order in 1/N C . The first steps towards 
establishing such an equivalence are just recently being taken P412]]. 

In either approach, a particularly fruitful physical process to examine has been meson- 
baryon scattering in the large- N c limit. The present paper furthers this study, taking as a 
tractable example of a multi-channel Lagrangian a variant of the linear cx-model. Before we 
specify the model, and our particular treatment of it, it is helpful to put the present work 
in historical context. 

A review of the relevant theoretical literature over the past decade reveals an interesting 
sociological phenomenon: there are two disjoint bodies of large-iV c papers devoted to two 
topologically distinct sets of diagrams, namely Compton-type versus exchange-type graphs, 
that contribute to the meson-baryon S'-matrix.[| Examples of Compton-type and exchange- 
type graphs are displayed in Fig. 1 and Fig. 2, respectively. Topologically, they differ in the 
following way: in the exchange-type graphs of Fig. 2, it is possible to trace a continuous line 
from the incoming to the outgoing meson without ever traversing a baryon line segment, 
whereas in the Compton-type graphs of Fig. 1 this cannot be done. 

Let us review, briefly, some of the salient points of physics that emerge from the study 
of each of these two classes of diagrams. 



So far as we are aware, the only attempt to date to treat these two classes of graphs in a unified 
manner can be found in Sec. 7 of Ref. Iiq|. 
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A. Compton-type graphs 



While presently the Compton-type graphs (Refs. [ I I ]l7j) are much less well understood 
than the exchange-type graphs discussed below, they nevertheless yield some interesting 
physics, as follows. Look at Figs, la and lb. Since each vertex scales like ^/W c (see Ref. ||), 
these graphs individually scale likeQ N c . However, we know from Witten's analysis of quark- 
gluon diagrams || that the total amplitude for meson-baryon scattering must scale like iV°, 
not N c . Therefore there must be leading-order cancellations between Figs, la and lb. Add 
to this observation another important piece of large- N c physics: the fact that for the case of 
two light flavors (which we focus on exclusively herein) the spectrum of stable baryons is a 
tower of states of equal spin and isospin [Qj: I = J = 1/2, 3/2, 5/2, ■ ■ • , N c /2, which are all 
degenerate in the large- N c limit (more precisely, in the limit J 2 /N c — > 0). We then demand 
leading-order cancellation between Figs, la and lb, for the reason described above, with 
the three baryon legs drawn from all possible baryon states in the I = J tower, consistent 
with triangle inequalities for isospin and angular momentum at each vertex. This exercise 
is carried out in Refs. |13| and | 14j . The upshot is a set of proportionality relations between 
the various coupling constants g n NN, 9-kNA, 9-kAA, and so forth up the I = J tower, relating 
each of these a priori independent couplings to a single underlying coupling constant, up 
to multiplication by Clebsch-Gordan coefficients. We call this set of relations for the pion- 
baryon couplings the "proportionality rule." Furthermore, Dashen and Manohar have shown 
that corrections to the proportionality rule do not occur at order 1/N C , as naively expected, 
but rather at order 1/N% This suggests that the proportionality rule should be relatively 
robust. Calculationally, it implies that, once the order N c contributions to the amplitude 
have cancelled, the surviving order A? pieces arise solely from the 1/N C corrections to the 
baryon propagator, and not from 1/A C corrections at the vertices, as one might have thought. 

Numerically, the proportionality rule for the pion-baryon couplings works well. Not only 
does the decay width of the A work out to within a few GeV of its measured value when 
g n NA is related, using this rule, to the experimental value of g n NN HEHi but furthermore, 
with the same input parameters, the widths of the "large- N c artifacts," i.e. the baryons 
with I — J > 5/2, are so large that they cannot be considered "particles" at all, and 
as such, pose no problem for phenomenology ||10||. This latter observation removes what 



2 The baryon propagator is approximated by i[v ■ k + ie) -1 in the large- N c limit, where v is the 
baryon's 4-velocity, k is the momentum imparted by the incoming meson (assuming the incoming 
baryon to be on shell), and it is also understood that one throws away the two small components 
of the Dirac 4-spinor. We focus on the kinematic regime k ~ N® so that the baryon propagators 
do not affect the N c counting. 
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has been, till recently, one of the chief objections to the entire large- iV c program. Another 
success of large N c is that the group-theoretic predictions of the old SU(2N F ) symmetry are 
recaptured P|6|, |l4l , |IE| , without the need to appeal to the construct of the nonrelativistic, 
weakly interacting constituent quark model. 

A further refinement was made recently by Jenkins ||17|| , who examined the one- loop 
chiral corrections to the masses Mj of the I = J baryons, and deduced the consistency 
relations 

Mj = Mo + + 0{N- 2 ) (1) 

where Mq and X are constants of order iV c that can be fixed, for example, by pegging M1/2 
and M3/2, respectively, to the experimental nucleon and A masses. 

While the large- N c results of Refs. [|l4l -|T7|1 are derived using effective Lagrangians of 
mesons and explicit baryons, the physics of the Compton-type graphs can also be accessed 
using the skyrmion approach . The parallelism between the two approaches is manifest 

in expressions such as Eq. ([J). In the language of the two-flavor Skyrme model, M and X 
are interpreted as the mass and moment of inertia of the soliton, respectively flSUl8|. It is 
reassuring that the expression ([]]) can also be gotten directly from looking at quark diagrams 
in large- N c QCD closing the circle. 



B. Exchange-type graphs 

Next we turn to the physics of the exchange- type graphs (Refs. |p||i9|-p7ll), which is the 
primary focus of this paper. Examples are shown in Fig. 2. These graphs likewise contribute 
to the scattering amplitude starting at order N®. Although the summation of all such graphs 
would appear to be an impossible task, it can actually be carried out in a straightforward 
manner — so long as one contents oneself with the leading-order answer in the 1/N C expansion 
0. As will be reviewed in detail in the Sections to follow, the key idea is to rewrite these 
multiloop graphs as trees, exploiting the large- N c approximation. Tree graphs have the great 
advantage over loops that they can all be summed by solving classical equations of motion. [] 



3 To remind the reader || that he or she already knows a situation where "loops" become "trees," 
recall the ancient problem of electron-proton scattering in the low-energy regime where the proton 
mass is much greater than all other scales in the problem. On the one hand, these are evidently 
multiloop interactions, in which the proton and electron lines exchange a large number of photons 
in all possible tanglings. On the other hand, we know that the physics is accurately described 
by classical equations: first the proton generates a classical Coulomb field, and then the electron 
propagates linearly through this non-trivial background (Rutherford scattering). These two dis- 
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It is this summability property which justifies our earlier statement that the exchange-type 
graphs are much better understood than the Compton-type graphs. 

While the analysis of this paper will be carried out using explicit baryon fields, the 
set of classical equations that emerges is, once again, highly reminiscent of the skyrmion 
approach, in which the corresponding classical equations describe a pion propagating through 
the background field generated by the skyrmion itself [|ITS|-|2"1|| . In particular, the group- 
theoretic relations familiar from the Skyrme model carry over intact to models such as the 
present one with explicit baryons. These include non-trivial, and experimentally reasonably 
well satisfied, relations in which isospin-3/2 ttN scattering amplitudes are expressed as 
linear combinations of the isospin-1/2 amplitudes fl9j22|| . Similar relations hold for kaon- 
nucleon scattering p3| , and for 7r N — ► pN |24] , and in fact for all quasielastic meson-baryon 
scattering processes. 



If, extending Donohue's original suggestion ||25|| , one crosses these relations among scat- 
tering amplitudes from the s-channel to the t-channel (e.g., NN — > mesons), they can be 
re-expressed concisely as two large- iV c selection rules [^,0. First, there is the very same 
"proportionality rule" discussed earlier, in the context of the Compton-type graphs. How- 



ever, the derivation given in Ref. p6| makes clear that the proportionality rule is completely 
independent of the chiral limit, and furthermore, that it applies not only to the pion-baryon 
couplings but equally to the baryon couplings of all bosons. Beyond the width calculations 
noted above P,pH|, the phenomenological validity of the proportionality rule is put to the 
test in Fig. 7 of Ref. P2| , in which the appropriate linear combinations of the experimental 
ttN — * tcN and ttN — > ttA scattering amplitudes are compared. 

In addition, a second large- N c selection rule emerges, the ll I t = Jt rule" [p6| , ^7|l . This rule 
states that the isospin of the emitted/absorbed meson must equal its total (spin + orbital) 
angular momentum, measured in the rest frame of the large- N c baryon. Concrete examples 
of meson-nucleon couplings that satisfy the I t = Jt rule include the pseudovector coupling 
of the pion, the tensor coupling of the p, and the vector coupling of the u p7) : 



(g M /2M N )d^^N^rN , g^ s d,p v ■ Na^rN , g™^ ■ N^N , (2) 

each of which must be augmented by couplings to the entire tower of I = J baryons as 
required by the proportionality rule. Since these couplings obey the I t = Jt rule, the three 
coupling constants are nonvanishing at leading order in the large- N c expansion: 



parate pictures are reconciled by the fact that the loop graphs are really trees, by exactly the 
same manipulations described in Sec. II below. The insight of Ref. || is that this same mechanism 
(modulo nonlinearities due to the fact that bosons, unlike photons, are self-interacting) holds for 
the exchange of arbitrary bosons in the large-iVc limit, thanks to the proportionality rule as well 
as the It = Jt rule reviewed below. 
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9itNN tens „vec 

2M, 



N 



gT ~ € c ~ • ( 3 ) 



In contrast, the It = Jt rule forbids at leading order the other two canonical vector-meson 
interactions, the vector coupling of the p and the tensor coupling of the uj, 

g v p ec p,-N^TN and g^d^ ■ Na^N , (4) 

meaning that these coupling constants must be down by (at least) one power of 1/N C com- 
pared to Eq. (0): 

„vec „tens ^ /-\ 

g P ~ — • ( 5 ) 



The relative unimportance of the vector (tensor) coupling of the p (u>) has long been known to 
nuclear physicists who construct one-boson exchange models of the nucleon-nucleon potential 
|p8|-pT| . It is pleasing to see these phenomenological rules of thumb emerge as theorems in 
the large- N c limit. 

C. Two interesting unresolved questions 

In lieu of a Conclusions section, we close this expanded Introduction with two questions 
that are food for further thought. First, is the complete meson-baryon S'-matrix properly 
obtained by adding the Compton-type and exchange-type graphs together, or, as an alterna- 
tive prescription, might it be the case that either set of graphs by itself (assuming an infinite 
spectrum of mesons) contains the complete answer? This latter possibility is suggested by 
the observation that mesons and baryons are composite particles made up from quarks and 
gluons. Since at the quark-gluon level there is no longer a topological distinction between 
the graphs of Fig. 1 and Fig. 2, one must be especially careful to avoid double counting, and 
this might conceivably preclude adding the graphs of Fig. 1 and Fig. 2 together in a naive 
way.Q 

Second, the exchange-graph formalism of Ref. || applies not only to the meson-baryon 
system which we focus on here, but equally to the baryon-baryon, baryon-antibaryon, 
baryon-baryon-baryon, and in general to all n-baryon, m-antibaryon interactions (Fig. 3). 
Of course, there are no analogs of Compton-type graphs for these multi-baryon systems. 
It follows that the exchange-graph formalism of Ref. || gives — in principle — the complete 



4 For the resolution of similar issues in atomic physics, namely the avoidance of double-counting 
when bound states are involved, see Ref. and references therein. 
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answer for these cases, to leading order in 1/N C . By this we specifically mean the following: 
given an effective hadron Lagrangian whose meson-baryon couplings properly embody the 
I t = Jt and proportionality rules, the complete set of Feynman diagrams of the sort exhib- 
ited in Fig. 3 can be summed to leading order in 1/N C . It would be an interesting exercise to 
carry out this program, starting from a well-motivated effective Lagrangian, and to compare 
the results to the popular Bonn [3~0| and Paris |H| potential models (which are derived from 
just the ladder diagrams with at most one crossing) as well as to the recent work of Weinberg 
and others that relies exclusively on chiral perturbation theory p4|. 



D. Outline of paper 



The remainder of this paper is organized as follows. In Sections II and III we review 
the exchange-graph formalism of Ref. |§ and apply it to two warm-up problems, a "cr-only" 
and a "7r-only" model, respectively. Sections IV- VI explore in detail the meson-baryon 
S-matrix in a richer model comprising both pions and a mesons, a variation on the Gell- 
Mann-Levy cr-model [Q. Obviously, we do not take this model seriously as a realistic 
depiction of hadron physics. Rather, we aim only to illustrate how the formalism of Ref. || 
leads in a concrete way to a quantitative calculation of the exchange-graph contribution to 
the multi-channel meson-baryon S-matrix. With the present model solved, the scene is set 
for more ambitious, realistic calculations, necessarily incorporating vector mesons. We are 
also interested in comparing the large-A c effective Lagrangian approach that uses explicit 
baryons, with earlier large-A c results from the Skyrme model. We come to the conclusion 
that much of the detailed structure of the meson-baryon S"-matrix which hitherto has been 
uncovered only with skyrmion methods, can equally be described by models with explicit 
baryon fields. At the same time, both approaches share significant problems in the low 
partial waves, the complete resolution of which remains a major technical hurdle. 



II. FIRST WARM-UP PROBLEM: A cj-ONLY MODEL 

As a first calculation, let us consider a model with only a mesons and (non-strange) 
baryons [pjp5|. Because the a has I = J = 0, this toy model avoids the spin and isospin 
complications due to non-commutativity of Pauli matrices at the meson-nucleon vertices. It 
also avoids the complications of inelastic 2-body channels (e.g., nucleons cannot turn into 
A's). 

The Lagrangian to be solved in this Section is the large- N c version of: 
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£*N = \{d^f ~ V(a) + jV(z 7 • d - M N )N - gaNN , (6) 
where, for definiteness, the a self-interactions are described by the fourth-order potential 

V{a) = ~m CT V + \™ Z + ^ ■ ( 7 ) 

By the words "large- N c version of" we mean that the coupling of the a to the nucleon in 
Eq. (|6|) must, in principle, be augmented by analogous couplings to the entire I = J tower of 
large- N c baryons, starting with the A (/ = J = 3/2) and continuing through the state with 
I = J = N c /2. The relative strengths of these couplings is given by the proportionality rule 

. However, in this simple model, since the a carries the quantum numbers of the vacuum, 
it couples diagonally to this tower (as noted above). Therefore, so long as we restrict our 
attention to nucleon targets, we can safely drop these additional couplings to the higher 
baryon states and work with the simplified Lagrangian (|B]). 

In the large- A^ c limit the nucleon has mass of order N c and its degrees of freedom freeze 
out. This means that the nucleon kinetic energy term in Eq. @ can be dropped, and 
the Yukawa term has iViV replaced by a static source j(x). The formal derivation of this 
intuitive prescription was given in Ref. ||. For completeness, we review it here. Looking at 
Fig. 4, the product of the nucleon propagators (reading from bottom to top) is 

i i i 
■x- x ■ ■ • x 



j> + # x - M N + ie i> + jfj + # 2 - M N + ie * j> + ^ H + ty n - X ~M N + ie 

~ 1 X---X- : . (8) 

2 fcio + ie fcio H h fcn-i,o + 

In the above we have taken the large- N c (i.e., nonrelativistic) limit of the nucleon propagators 

i> + ^ - M N + ie 2 k + ie ' U 

assuming that the nucleon is in its rest frame. The prefactor (70 + l)/2 is the projector 
onto the large components of the Dirac 4-spinor. From now on we suppress it, with the 
understanding that we always throw away the small components. 

Our desired result is obtained by summing over the n! crossed ladders (Fig. 5), and using 
the interesting identity for distributions, 

n i i i 
2n5( k i0 ) x x • ■ ■ x 

2=1 permutations "'ilO 

+ ie k h0 + k i2 o + ie k h0 H h k in _ l0 + ie 

(«l.—>«n) 

= 2-K5(k 10 ) x 2TT5(k 20 ) x ■ ■ • x 2-K5(k n0 ) . (10) 

(To prove this, Fourier transform both sides of this identity in all n momenta.) Each of the 
n\ terms in this sum corresponds to a distinct crossing or ordering of the n rungs of the 
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ladder. The 5-function on the left-hand side of this equation reflects conservation of energy 
along the nucleon line in the large- N c limit: 

n n 

2Tr8(-p' + Po + Y,kia) ^ 27nJ(S>o)- (11) 

i=i i=i 

Recognizing 27r<5(A;o) as the 4- dimensional Fourier transform of <5 3 (x), we immediately 
understand the meaning of the simple factorized right-hand side of Eq. ([IT]) in terms of 
graphs. Simply put, the sum of the n! crossed ladders is equal to the single graph of Fig. 6, 
generated by the effective Lagrangian 

£ eff = ^(d»*) 2 -V{*)-aj{x) (12) 

where, as promised, the nucleon field has been frozen out in favor of the external c-number 
source 

j(x) = ^ 3 (x) . (13) 

The complete exchange-graph contribution to <jN scattering in the large- N c limit now 
emerges from a two-stage numerical program, which is most transparent in graphical terms. 
In the first stage, one defines a "classical" a field a c \ as the sum of all one-point trees 
(Fig. 7). The reason one considers only trees is that meson loops are suppressed by powers 
of l/N c @,H,|J. In the second stage, one considers a propagating a field (which we call the 
"quantum" field cr qu to distinguish it from a c \) interacting with an arbitrary number of a c \ 
insertions (Fig. 8). By inspection, this two-stage procedure is equivalent to summing all the 
tree graphs of the form shown in Fig. 6 (the loop graphs being sub leading in 1/N C ). As 
promised: the loops (Figs. 4-5) have turned into trees, exactly as in the old electron-proton 
problem invoked in Sec. I. 

This two-stage graphical procedure is easily translated into the language of differential 
equations. Solving for a c \ as per Fig. 7 is equivalent to solving the classical Euler-Lagrange 
equation for the effective Lagrangian of Eq. flT2]), namely, 

- W cl (x) + y'(a cl (x)) + j(x) = . (14) 

Note that a c \ is time-independent because the source j(x) has this property. Next, solving 
for the propagating field cr qu , as given by Fig. 8, is accomplished by noticing that at every 
vertex, there are exactly two cr qu legs, the rest being insertions of a c \, with the coupling 
constants read off from V(a). Therefore, the relevant equation of motion comes from the 
quadraticized Lagrangian 

£quad = ^Oqu^Oqu ~ ^qu^'Vclf*)) , (15) 
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which induces the linear time- dependent equation 

[d^ + U> cl (x))] a qu (x) = . (16) 

In short, we have outlined a two-stage numerical procedure, the first stage involving a 
non-linear time-independent equation for a "classical" meson field, the second involving a 
linear time- dependent equation for a "quantum" meson field in the classical background. 
This is reminiscent of the skyrmion approach to meson-baryon scattering Hl9|-|2^|. In the 
subsequent Sections, when pions are introduced, this correspondence will be sharpened by 
the emergence of a hedgehog structure to the classical pion field that is familiar from the 
Skyrme model []7]||. (The chief difference between the two approaches is, of course, that 
baryon number is carried by topology in the Skyrme model, and by smeared 5-function 
sources when the baryon fields are explicit.) 

The analog of the hedgehog Ansatz in the present model with I = a mesons alone is 
just ordinary spherical symmetry: 

<r cl (x) = G{r) . (17) 
The profile function G(r) is found by solving the non-linear radial differential equation 

G" + -G' - m 2 G - £ G 2 - ^ G 3 = j(x) (18) 

implied by Eqs. (^) and (fHI). Unfortunately, Eq. (0) suffers from ultraviolet problems 
when j is literally taken to be a ^-function as per Eq. ([13]). The source of these divergences 
(which are worse than in the original loop graphs, Figs. 4-5) can be traced to the nonrel- 
ativistic reduction of the propagator (||), which is only valid so long as the components of 
the exchanged meson momentum satisfy <C M^. (A similar breakdown of the large- N c 
approach is discussed in Sec. 8 of Ref. 0.) A simple cure is to smear out the source, say, as 
a Gaussian: 

j(x) — > 9 exp(-r74) , a N ~N° c . (19) 
(ajvV^J 

This approximation now renders Eq. ( |T8"D tractable, at the expense of introducing a "nucleon 
size" parameter into the problem. This new parameter provides an ultraviolet cutoff on 
the momentum allowed to flow into or out of the nucleon. We have checked that our 
numerical results are not overly sensitive to over a reasonable range of values. 

Equation ([18|) represents a two-boundary value problem which can be solved in an itera- 
tive fashion using a standard "shoot and match" Runge-Kutta integration procedure p6|.p| 



5 For details, see the Appendix. 
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Fig. 9 shows the profile function G(r) for the specific choice of parameters m CT = 600 MeV, 
k = 18.5, A = 214, g = 13.6 and a at = 0.5 fm. Note that G(r) looks very much like a 
Yukawa function, exp(— m a r)/r, except that it is finite at the origin (because of the smear- 
ing of the nucleon source term) and has small deviations in the 0.5 to 1.0 fm region due to 
the non-linear terms involving k and A. 

Given G(r), we then solve Eq. (^) for cr qu by means of a standard partial wave analysis. 
For angular momentum /, the radial scattering wave function Ui(r) = r<Ji(r) having energy 
us satisfies 



+ q - kG(t) G z (r) 



dr 2 2 r 2 



ui{r) = 0, q 2 = uj 2 -ml. (20) 



This is a Schrodinger-like linear differential equation that can also be solved by Runge-Kutta 
integration from the origin (where ui(r) must be regular, going like r l+1 ). The asymptotic 
form of ui(r) is then fit in the usual way to a linear combination of spherical Ricatti-Bessel 
functions, ji(qr) and ni(qr), yielding the phase shifts for aN elastic scattering. 

The potential in Eq. ([TBI) [or Eq. (|20|)1 has a short-range repulsive core (coming from the 
quartic term in the Lagrangian) and intermediate-range attraction (coming from the cubic 
term and the fact G(r) < 0). Consequently, as shown in Fig. 10, the S'-wave phase shift at 
low energies is positive because of the medium-range attraction, but it soon turns over and 
looks like the phase shift for a hard-core repulsive potential. At still higher energies (not 
shown), the phase shift returns to 0, since the short-range repulsive core is finite. The higher 
partial waves exhibit similar behavior, but offset to increasingly higher energies because of 
the angular momentum barrier. 



III. SECOND WARM-UP PROBLEM: A vr-ONLY MODEL 

As a second simplified example, we consider a model of gradient-coupled pions and I = J 
baryons. The Lagrangian we want to solve is the large- N c version of 

f ^ £1 ^ a PiV-^a- "^™2 —a^a ^ I _a„a\2 

C wN = -C M vr c^vr - ~m n n n - — (tt tt ) 

+]V(i7 • d - M N )N - g^n ■ N^^tN . (21) 

The reason for choosing pseudovector coupling rather than the pseudoscalar coupling, —g' n TT- 
Nij 5 TN, is that it is more amenable to a large- N c treatment, for the following reason. The 
matrix 7 5 is purely off-diagonal, connecting the large to the small components of the Dirac 
4-spinor. This means that taking the nonrelativistic limit of the baryons is a singular 
operation when the pion is not soft. In contrast, 7 5 7 M does connect the large components 
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to themselves for fi = 1, 2, 3 so that, with a pseudovector coupling, we can follow the simple 
leading-order large- N c prescription given earlier of just throwing away the small components 
(including, inter alia, the 7 5 7° contribution; remember that the 1/N C expansion breaks 
up Lorentz invariants). Of course, in a different limit from large N c , namely the soft- 
pion limit in which the pion is emitted from the on-shell nucleon at approximately zero 
4-momentum, the pseudovector and pseudoscalar couplings are indistinguishable, provided 
one takes g w = g' n /2M N . 

The meaning of the words "large- N c version of" preceding the Lagrangian (|2~1~D is that, 
as in the aN model of the previous section, the coupling of the pion to the nucleon must 
be supplemented by analogous couplings to all the other members of the tower of I = J 
baryons, and likewise for the nucleon kinetic term. In the previous case this was an irrelevant 
complication: because the zero-isospin a cannot induce transitions between states in this 
tower, the problem diagonalizes. In contrast, pions can and do change nucleons to A's, etc. 

The most convenient way to implement the gradient coupling of the pion to the I = J 
tower of baryons is to change baryon basis to the so-called collective coordinate basis \A) 
familiar from the Skyrme model, with A an SU(2) element ||. These basis elements are 
defined by [g,[18| . 



\A) = £ (2R + 1)^ £ H^D^liM^) A ) . (22) 

fl=l/2, 3/2,- i z ,s z =-R,-,R 



normalizing the volume of SU(2) to unity. On the right-hand side of this equation, the 
baryons are given in the usual spin-isospin basis, e.g., a neutron of spin up and a A + of spin 



projection —3/2 would be denoted as i/^i/^) an d 1/2-3/2 )> respectively. In the collective 



-1/2,1/2 



3/2 



1/2,-3/2/' 



coordinate language, the correct pion-baryon coupling reads 

- 30, £ d b n a I dAD%(A) \A) (A\ . (23) 

JSU 2) 



o,6=-l,0,l 



This coupling was first written down on general grounds (without reference to soliton 
physics) by Adkins, Nappi and Witten ||, and is necessary for the consistency of the 
Compton-type graphs with the overall N® scaling of the pion-baryon scattering amplitude 
in the large- iV c limit, as reviewed earlier [T^JT^]. It has also recently been established using 
collective coordinate quantization of the skyrmion |10|| . Despite our convenient adoption 
of Skyrme-model notation, we emphasize that the coupling (^) is, in the present context, 
understood to be built from explicit baryon field operators, and not solitons. 

Let us verify explicitly that Eq. ( p3|) is indeed the correct large- iV c pseudovector coupling 
of the pion to the baryon tower. In particular, Eq. ( p3l) has the following four desirable 
properties: 
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1. It is invariant under isospin and angular momentum; 

2. It contains the pion-nucleon interaction shown in Eq. (Ell); 

3. It correctly implements the "proportionality rule" governing couplings to the higher 
states in the I = J tower. 

4. It accurately predicts the width of the A, and furthermore, gives widths so large for 
the large- A/" c artifacts of the model (the baryons with / = J > |) that these pose no 
phenomenological problems for the large- N c approach. 

We deal with each of these assertions in turn: 
1. The state \A) transforms as 



\A) ^ \UjA) and \A) 



so that 



SU(2) 



dAD$(A)\A)(A\ 



Jsui2)dAD%(A) 



AU] 



U T AU\)(U T AU 



(24) 



Ism^dAD^UjAUj) \A) (A\ 



)SU{2) 

D$(V})D&(rfj) J sm dAD%(A) \A) (A\ 



(25) 



Here we have used the group invariance of the SU{2) measure, d{U\AUj) = dA, and the 
reality property of the D {1) matrices. Similarly, 



d b 7T a 



(26) 



Combining these last two equations and using the composition property of the Wigner D 
matrices, we confirm that the coupling of Eq. (|23|) is invariant under isospin and angular 
momentum rotations. 

2. Using Eq. (0), we rewrite the coupling of Eq. (p3|) as 

R' " 



R,i z ,s z FCi'^s'z 



a.b 



X 



SU(2) 



) i? '"^[(2 J R+l)(2i? / + l)] 1/2 
dAD^{A^%AA^)D^,M ] ) 



R 



-3g n J2d b n a E E (-) R+R> (Rli z a\R' i' z )(R' ls' z b\Rs z ) 

a,b R,i z ,s z R',i' z ,s' z 

R' \ I R 



x 



1? s 7 



(27) 



using standard D-matrix integration tricks. We now pick out the terms with R = R' = 1/2 
in this expression in order to study specifically the pion coupling to iViV. Isospin and angular 
momentum invariance can be made more manifest by rewriting this subset of terms as 
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S^E E E 



1/2 



a, 6 i z 



'1/2 



(28) 



which we recognize as the nonrelativistic (or, equivalently, in the present context, large- N c ) 
limit of the gradient coupling —g^d^n ■ N^j^tN. 

3. A careful reading of Ref. [^6] reveals this criterion will be automatically satisfied due 
to the diagonality of the pion-baryon coupling, Eq. (23), in the collective coordinate A. It 
is instructive nevertheless to see how this comes about explicitly The baryon-antibaryon 
Hilbert-space operator in Eq. (|27|) can be written in terms of states with good t-channel 
(exchange-channel) quantum numbers as follows: 



R! 



R 

i z s z 



E E n R+ ' z n R,+s Hhi tz \R'R^^ z )(RR's z ,-s:\j t j tz ) 

hjtz Jt,Jtz 



X 



It 5 RR' \ I Jt \ RR' 



Ltz 



t 1 
Ji 



tz 



(29) 



where the phases in the above are the usual cost of turning bras into kets in SU(2) 



\jm) 



-m 



Plugging Eq. (|29| ) into Eq. ( p?| ) and using Clebsch-Gordan orthog- 



onality gives for the pion-baryon coupling: 

h. X- (-\ R + R '[(2R +1) (2R' + 1)] 1/2 



.'/•- E 

Itz,Jtz 



E 

R,R' 



l;RR n 

hz , 



Jt 



-l;RR' 

Jtz 



(30) 



This equation correctly embodies two large- N c selection rules: the fact that the exchanged 
angular momentum J t is equated to the isospin I t — 1 of the pion is a specific example 
of the more general I t = Jt rule [p6,27f|, whereas the square- root proportionality factors 
relating the pion's couplings to the various baryon states in the I = J tower illustrate the 
proportionality rule [p6fl . 

4. The coupling (|23|) can be used to calculate the decay width of a baryon with 
spin/isospin J to the next-lower state J — 1 via the emission of a single pion. For the 
case A — > Nn one calculates Ta = 114 GeV as against a measured width of 120±5 MeV 
Pleasingly, for the higher states, 



1C 



/ = J > |, the widths turn out to be so large 
that these large- N c artifacts cannot be said to exist as particles, and therefore, pose no phe- 
nomeno logical problem for the large- N c program. One finds Ts pa 800 MeV, Tz pa 2600 MeV, 



6400 MeV, and so forth H0f 



As before, we seek to sum the set of exchange-type graphs of the form shown in Fig. 2. 
However, a priori, the situation is not so simple as in the a-only model of Sec. II. Look 
again at the interesting identity ([II]) for distributions, which is the key to turning loops into 
trees. The n\ terms on the left-hand side correspond to the n\ distinct "tanglings" in which 
the exchanged a lines are attached in a different order to the baryon line. Because the a 
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carries no spin or isospin, each tangling enters with the same relative group-theoretic weight 
in Eq. (TO), and the identity goes through as written (so, too, for photon exchange). In 
contrast, tt, p and uj mesons, etc., carry non-trivial isospin and/or spin, and the n\ tanglings 
would not be expected to occur with the same group-theoretic factors. (Pauli spin/isospin 
matrices do not commute.) Specifically, one expects a different product of n spin and n 
isospin Clebsch-Gordan factors weighting each term on the left-hand side of Eq. (PTOQ , and 
destroying the identity. Nevertheless, acting together, the I t = Jt and proportionality rules 
assure that, to leading order in 1/N C , these n\ group-theoretic factors are indeed equal, once 
the intermediate baryon legs are summed over all allowed I = J states. Therefore, the 
identity fllPf) , derived for a (or photon) exchange, applies as well to the exchange of these 
non-trivial mesons. This theorem is proved in Ref. @, using elementary properties of 6j 
symbols. However, there is an easier way to see this, which is to work directly in the \A) 
basis. So, look again at Fig. 5, and understand the baryon line to mean, not a nucleon or 
a A or any specific member of the I = J tower (which can change identity at each pion 
interaction vertex), but rather a baryon state \A) sharp in the SU(2) collective coordinate 
A, which is preserved at each vertex, due to the diagonality in A of the coupling (p3|) . Initial 
and final nucleon, A, etc., states can be projected out at the very end of the calculation 
using standard group-theoretic techniques borrowed from the Skyrme model [i.e., inverting 
Eq. (p2|)]. At earlier stages, however, we can use the full machinery of Sec. II to turn loops 
into trees with impunity. 

Therefore, once again, the graphs of Fig. 5 can be summed following a two-stage program. 
In the first stage, one solves a static non-linear equation for 7r c \(A) (noting that the classical 
pion field depends on the SU(2) collective coordinate A). Isospin covariance trivially relates 
this quantity to n c i(A = 1), henceforth called just 7T c i. Using D^(A = 1) = 5 a b, one obtains 
the Euler-Lagrange equation 



by assuming a hedgehog Ansatz for the classical pion field (anticipating the resemblance to 
the Skyrme model): 



subject to the boundary conditions that F(r) be regular near r = and bounded as r — > oo, 





(32) 



Equation (|3~1~D then becomes an ODE for the classical pion profile F(r): 




(33) 



15 



F(r) = Br + C(r 3 ) near r = 0; 

F(r) — > Cexp(— m 7r r)/r as r — > oo . (34) 

5 and C are scale parameters that are initially unknown to us but are fixed implicitly by 
the non-linearity of Eq. (|33|). This is another two-boundary- value problem, which can be 
numerically solved as before (see Fig. 11, and Appendix A). 

In the second stage, one solves the linearized time-dependent equation for 7r qu propagating 
in the background of tt c \(A). Again, isospin invariance trivially relates this process to the 
propagation of 7T qu in the background of 7r c i(v4 = 1), the latter quantity being given by 
Eqs. (|3^) and (^). Initial and final nucleons or A's are then projected from the hedgehog 
by inverting Eq. fl22|) , using the orthogonality over SU(2) of Wigner D-matrices. Finally, 
the initial and final pion-baryon systems are combined into states of good total isospin and 
angular momentum in the usual fashion to give the partial-wave S'-matrix for ttN — > ttN, 
ttN — > 7rA, etc. 

Fortunately, this cumbersome (if straightforward) sequence of group-theoretic steps can 
be circumvented, once one realizes that they are identical to the procedure followed in the 
Skyrme model [|I~9|-|22||. Rather than "reinventing the wheel" one can therefore carry over 



intact the machinery of Refs. p9|-p^1 of i^-spin decomposition and 6j symbols. We postpone 
the explicit review of this formalism to Sec. IVTl in which we complete the analysis of the 
richer model containing both pions and a mesons. 

Unfortunately, the pion-only model discussed in this Section is inherently uninteresting 
phenomenologically. Because of G-parity, the pion-pion interactions can only come from 
even powers of n(x), which means that the potentials entering into the coupled Schrodinger- 
like scattering equations are strictly repulsive. [They are proportional to XF 2 (r).] As a 
result, all nN phase shifts exhibit repulsive behavior (i.e., clockwise motion in the Argand 
plots with increasing energy). Thus there is no possibility for tcN resonances in such a 
model. We need something like the a meson to provide a range of attraction between 7r's 
and JV's. 



IV. DEFINING THE cr-vr MODEL 

In view of the two models discussed in the two previous Sections, one might have some 
hope that a model combining a and 7r mesons would provide a more promising (if still 
crude) description of pion-nucleon interactions. In this model the a meson will be taken as 
an "elementary" field, along with the three ir fields. Indeed, in the large-iV c limit, the a, if 
such a state exists, is necessarily a stable particle, as the decay amplitude to two pions is 
suppressed by l/y/N c . 
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For guidance in constructing our \arge-N c model of pions and a mesons, and selecting 
reasonable values of the coupling constants, we recall the linear cx-model of Gell-Mann and 
Levy: |2| 



£ = 1<V«9V + i^Tf • d^n - ^{a' 2 + tt • tt - a 2 ) 2 + aa' 

-ga'NN - gn ■ Ni^rN + Ni-y ■ dN . (35) 



In this well-known model, the nucleon and a get their masses through dynamical symmetry 
breaking, the a vacuum expectation value v being g~ x M^, and chiral symmetry emerges in 
the limit a — > 0. It is convenient to redefine the a field by subtracting the VEV, 

cr'(x) = v + a(x) . (36) 

By substituting for o' and expanding, the four coupling constants {g, A, a, a} can be traded 
for the more physical set of parameters, {g, Mjv, m^, m^}, using 

, S 2 , 2 2l m, 2 M N 2 M% (ml - 3m|) 

= _ a = ^- 0= ^w^y (37) 

In this paper we will take 

g = 13.6 , M N = 5.0 far 1 , m w = 0.7 far 1 , and m a = 5.0 far 1 . (38) 

This choice for the nucleon mass roughly averages the actual iV and A masses, while the a 
meson here could be identified with the /o(975) meson for concreteness. The value of g is the 
measured pion- nucleon pseudoscalar coupling constant. With these values, the non-linear 
self- interaction strength has a large value, A ~ 91. 

For a large- N c treatment, the Gell-Mann-Levy model needs to be modified in the follow- 
ing two ways. First, as discussed in Sec. [TTJ, the pseudoscalar nN coupling is inappropriate, 
and should be replaced by pseudovector coupling as in Eq. (^l|), with g n = g/{2M N ) = 1.42 
fm. Unfortunately, with this replacement chiral symmetry is lost, even for a = 0. However, 
as stated in the introduction, our purpose in this paper is to explore the large- N c approach 
in a multi-channel model, not to present a fully realistic effective Lagrangian of the low-lying 
hadrons, which would require not only approximate chiral symmetry but also the incorpo- 
ration of vector mesons. (To look at the bright side, the fact that we are sacrificing chiral 
symmetry re-emphasizes the point that our large- iV c techniques have nothing to do with the 
chiral limit.) 

Second, the meson couplings to the nucleon must be augmented by suitable couplings to 
the entire I = J baryon tower (and likewise for the nucleon kinetic energy). The prescription 
for doing so is Eq. ( |2"3"D for the pion. It is easy to check that the analogous prescription for 
the a is given simply by 
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— gaNN 



-go 



SU{2) 



dA \A) (A\ 



(39) 



As previously, we solve for the classical meson fields, for the reference choice of collective 
coordinate A = 1, by means of a hedgehog ansatz: 



<i(x) 



r a F(r) 



G(r) 



(40) 



Smearing out the 5-function baryon source as in Eq. (|19"D, we find coupled non-linear Euler- 
Lagrange ODE's for F and G: 

2 



dr 2 r dr 



+ mt) F-X 



F 6 + FG + 2vFG 



3gr 



M N a%7r^ 



exp (^—r 2 /t 



N 



*G+-±G-mlG 
dr r dr 

9 



A 



G 3 + F 2 G + 3vG 2 
exp (-r 2 /a 2 N ) + XvF 2 



(41a) 



(41b) 



(a N ^) 3 

We will generally set the nucleon size parameter a at = 0.52 fm, but we will also consider 
the dependence of our results on in Sec. VI(C) below. The boundary conditions are that 
F and G must be regular at the origin and exponentially decaying (rather than growing) at 
infinity. The classical pion profile F(r) falls off like exp(— m n r)/r at large distances. On the 
other hand, G(r) falls off not like exp(— m cr r)/r as one might naively expect, but rather like 
exp(— 2m n r)/r 2 due to the F 2 source term on the right-hand side of Eq. (|41b|) , and the fact 
that 2m n < m a . Details of our numerical "shoot and match" procedure for solving Eq. ( |4"1"D 
can be found in Appendix A. 

The solution for F(r) and G(r) is shown in Fig. 12. Note that G(r) is negative with 
respect to F(r) and v. It is this relative sign that leads to the attractive nN interaction 
found in this model. 



V. PION-HEDGEHOG SCATTERING 



Having solved for the classical pion and a fields, we turn to the small-fluctuations problem 
of meson-baryon scattering. As in the Skyrme model p9|-p2|, one first solves for meson- 



hedgehog scattering, and subsequently one folds in some group theory (6j symbols) to obtain 
meson- nucleon scattering. The meson-hedgehog S- matrix is the topic of this Section, while 
the meson-nucleon S'-matrix is the subject of Section VI to follow. 

We return to the a-n Lagrangian, Eq. fl3"o] ) as modified subsequently in the text in the 
manner suggested by large- N c . Consider fluctuations of the meson fields about their classical 
solutions, 



18 



IT (X) 



r a F(r) + < u (x) , a(x) -> G(r) + a, 



qu 



(42) 



Since F and G satisfy the Euler-Lagrange equations, terms linear in the fluctuating fields 
vanish. The quadratic terms then lead to linear equations of motion for tt^ u (x) and cr qu (x). 
Higher-order nonlinearities in the meson fields are subleading in 1/N C , as previously noted. 

We will work out the partial-wave scattering amplitudes factoring out a uniform time- 
dependence exp(— iujt) from all the fluctuating fields. For the a this involves the usual 
expansion in spherical harmonics, 



a qu ^,Xj 



Yl <l>KK z {u,r) Y KKz (±) 



(43) 



K,K Z 



For the pions the decomposition is slightly more complicated ||19|--f22|] . The conserved quan- 
tum numbers are not isospin and total angular momentum but the so-called "grand spin," 
K = I + J. Since pions are spinless, J is just L, the orbital angular momentum. Thus the 
appropriate partial wave analysis for pions involves an expansion in terms of vector spherical 
harmonics, 



7T qu (^,x)= ^KK z L{uJ,r)y K L Kz 



(44) 



K,K Z ,L 

where L runs over values K — 1, K, and K + 1. 

For each value of K, the equations for the four radial wavefunctions 4>k, iPk,k, and 
ipK,K±i might be expected to form a 4 x 4 coupled system, [] but parity uncouples ipK,K from 
the other three. It obeys 

K(K+l) 



d 2 , 2d, 

-T^WK,K H r^K,K + 

dr z r dr 



K(r) 



^K,K = 



where 



uo 2 - ml and V v (r) = \[F 2 (r) + G{r){2v + G{r))\ . 



(45) 



(46) 



The remaining 3x3 coupled system of equations^ is best expressed in matrix form. 
Assembling ipK,K±i and <$k into the column vector 



VK,K-i{r) 



fl>K,K+i(r) 
<f>K(r) 



(47) 



6 From now on we drop the K z label on cj) and tp since the ensuing equations are independent of 



7 For the special case K = this is a 2x2 as ipo _i does not exist. 
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wc 



findQ 



^ K + ~* K + [Qjc -Vk]-*k = 0. (48) 



Here Qk is the diagonal matrix 

f 9 (K-1)K o (K + 1YK + 2) 9 K(K + 1)\ , . 

= dia S [£ ~ r 2 — . <fi " " £ ^ - r 2 J j , ( 4 9) 

and Vx is the symmetric potential energy matrix 

V u = V n (r) + 2XF\r) (j^-j) (50a) 



( . . K(K + 1)\ 
V 12 = -2AF 2 (r) \ V 2R+1 1 (50b) 

V 13 = 2AF(r)(^ + G(r))(^ T ) 1/2 (50c) 
V 22 = ^W + 2AF 2 (r)(^ r ) (50d) 

V 23 = -2XF(r)(v + G(r))(^±^j 1/2 (50e) 
V 33 = V a (r) , (50f) 

where we have defined 

g 2 =CJ 2_ m 2 and V a (r) = X[F 2 (r) + 3G(r)(2f + G(r))] . (51) 

Note that q% can be positive or negative, depending on whether the energy u is above or 
below the a threshold. 

The "diagonal" potentials V n and V„ are plotted in Fig. 13. They are repulsive at short 
distances and attractive at intermediate range. The factor of three in the definition of V a 
makes it about three times more repulsive and attractive than V n . Note that the vertical 
scale is in inverse fermis; these are potential wells of depths about 6 and 2 GeV, respectively, 
which means there is substantial attraction in both the aN and ttN systems. Also shown 
in Fig. 13 are the off-diagonal transition potentials V 12 and V 13 (but without iT-dependent 
factors) which are comparable in size to the diagonal potentials. 

Numerically, the uncoupled equation ( f45|) is readily solved using the Runge-Kutta tech- 
nique employed in Sees. II and III. This method also works for the coupled equations (HBT) , 



In so doing we are greatly assisted by the vector spherical harmonic identities, Eq. (10), in 



Ref. [21 1 . Note a typo there: K in the numerator of the square-root in the middle line of Eq. (10) 



shouldinstead be K + 1. 
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but only above the cx-threshold, 10 > m^. The problem below threshold is to ensure that the 
o wavefunction remains exponentially decaying, 



(r) — > exp(— Kr)jr , k = (m^. — u 2 ) 1 ^ 2 



(52) 



In our experience, numerical noise in the Runge-Kutta integration invariably induces ex- 
ponential blow-up: <pK{ r ) exp(+/cr)/r. We emphasize that even below threshold the a 
cannot be neglected as it causes substantial attraction in the nN channel. (Recall that the 
"box diagram" for ttN — > o~N — > ttN, Fig. 2c, is attractive.) 

A numerically more robust approach that works both above and below the a-threshold 
is to convert Eq. (f48|) into a set of coupled Fredholm integral equations of the second kind, 



*g>(r) 



Jk\t)+ I G K (r,r')\/ K (r')^(r')dr' 



(53) 



where the index i labels the linearly independent choices of inhomogenous driving terms. 
Above the a threshold, i runs over 1,2,3 and the inhomogeneous terms are 

( - \ / n \ 



J {1) (r) 



\ 








, J {2 \r) 



J 





j K +i(q«r) 




, J {3) (r) 



J 








(54) 



Below threshold, only the first two of these should be kept. The multi-channel Green's 
function Gk is the diagonal matrix 



Gn(r, r') 
G22O", r') 
G33O", r') 
G33O", r') 



— jK-l(?7r^<)^K-l(g7r^>) , 

—jK+i^r^nK+^ry) , 
2 A 

— £if(/«r<)A;if (ftr>) , below threshold , 

— jK(qa^<)nK(qar>) , above threshold . 
qa 



(55) 



where ji, ni are spherical Ricatti-Bessel functions ]38|] and %i, k t are modified spherical 
Ricatti-Bessel functions P5| , regular at the origin and exponentially decaying, respectively. 
By design, the multi-channel Green's function assures regularity of the wave functions at 
the origin and the asymptotic exponential fall-off of the 4>k below the a threshold. Note 
that G33 is continuous through the threshold. 

The S-matrix for the uncoupled pion scattering, Eq. (145|) , will be denoted here as the 
single-subscript quantity Sk, where the orbital angular momentum quantum numbers L = 
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V = K are suppressed. It is derived from the asymptotic analysis of the wavefunction in 
the usual way. The corresponding phase-shift 5k, defined as 



S K = e 2 ^ , (56) 

is plotted against pion momentum k in Fig. 14 for K < 5. For each K the corresponding 
phase shift is attractive, if numerically small apart from the case K — 1, and comparatively 
much less significant than in the Skyrme model (cf. Fig. 1, Ref. [^T[]). As always in scattering 
problems, the centrifugal barrier term in the scattering equations delays the onset of the rise 
in the phase-shift for the higher- L partial waves. 

The coupled-channels 3x3 (above threshold) or 2 x 2 (below threshold) part of the 
S'-matrix will be denoted Sfj, i, j = 1, 2 and/or 3, according to L — K — 1, K + 1, and/or 
K. It is obtained as follows. First, the K^-matrix is formed according to 

K£ = / drh{q 3 r)[\/{r)^\r)] 3 , (57) 

where L = K — l,K + 1 and/or K for j = 1,2 and/or 3, respectively, and also q\ — q 2 = q n , 
<?3 = Qcr- From the K^-matrix, the S- matrix is formed in the usual way, 

= feM) V2 [(l - iK*)(l + zK K )~% , (58) 

where for an explanation of the square-root flux factors (needed only for multichannel scat- 
tering) we refer the reader to Ref. [iO]. Time-reversal invariance implies S K = (S K ) T , 



which we have found to be a stringent check on our numerics. We will parametrize as 
r}f- exp 2i5^ subject to this symmetry property as well as to unitarity. 

The phase-shifts corresponding to the specific S'-matrix elements and are plotted 
in Figs. 15-16. Recall that, with our notation, these are the S'-matrix elements that describe 
pion-baryon scattering (no 'in' or 'out' er's, only intermediate cr's) in which the orbital 
angular momentum quantum number is preserved (L = L'), as opposed to changing up or 
down by two units (as it can for ttN — > irA). The bulk of the attraction in the present 
model, due primarily to the intermediate a-meson states, shows up in the phase-shifts of 
Fig. 16, with L = V = K + 1. Here one sees resonances (phase shifts rising rapidly through 
90 degrees) in each partial wave. In the L — L' — K — 1 partial waves (Fig. 15), one also 
sees attraction, although not so strong as to produce resonances. The surprise here is in 
the channel K=1,L = L' = 0, which reveals the presence of a bound state (Levinson's 
theorem). Once one folds in the appropriate group theory in the following Section to project 
the hedgehog onto physical baryons, the existence of such a bound state manifests itself as a 
parity conjugate to the nucleon. This feature is, unfortunately, not found in Nature, nor in 
the Skyrme model, and is an unwanted, unphysical artifact of the present strongly-coupled 
<7-7r model. 
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On the other hand, an improvement over the Skyrme model is the fact that all these 
phase-shifts (as well as those not plotted) eventually return to zero for sufficiently high 
energies. In contrast, in the Skyrme model they apparently grow without bound, eventually 
violating the unitarity constraints of quantum field theory, although admittedly at energies 
where several key approximations made in Refs. fll9|-|2"2|j , such as the neglect of skyrmion 
recoil, are clearly unwarranted. 

Another point to note about Fig. 15 are the cusp effects due to the opening of the a 
threshold at 10 = 5 fm~ . This is most apparent in the K — 1,L — phase shift, but the 
effect is present in the higher partial waves as well. 



VI. ttN elastic scattering 



A. Group-theoretics for meson-nucleon scattering 

In the previous Section we derived an ^-matrix for the scattering of pions and cr's off 
hedgehogs. The scattering information is encoded in partial-wave amplitudes we called Sk 
and Sfj where S K is a 2 x 2 matrix below the a threshold and a 3 x 3 matrix above it (except 
when K = in which case S K is 1 x 1 or 2 x 2, respectively). The integer index K labels 
the vectorial sum of the incoming or outgoing meson's isospin and angular momentum. K is 
conserved when the meson scatters off an object with hedgehog symmetry, in the same way 
that orbital angular momentum L is conserved in scattering from a spherically symmetric 
potential. 

Of course, what we are really interested is scattering, not from a hedgehog, but rather 
from a nucleon or A. The relationship between the two problems, "physical scattering" versus 
"hedgehog scattering," is contained in the following group-theoretic expression [|P9"|--f2~3|] : 



SW^Wto» = s kll'(uj) ■ (-) R '~ R [(2R + l)(2R f + 1)] 1/2 (2K + 1) 

K 

{K I tot Jtot\ [ K I tot Jtot I 
(59) 
R L I J [R' V V J 

Here u is the meson energy in the baryon rest frame (baryon recoil being subleading in 
1/N C ), L (L f ) is the initial (final) orbital angular momentum, / (/') is the isospin of the 
incoming (outgoing) meson, and R (R') is the spin/isospin of the initial (final) I = J baryon 
(e.g., R = 1/2 for a nucleon, R = 3/2 for a A, etc.). For physical scattering, K is no 
longer conserved; it is just a dummy of summation, constrained by the triangle inequalities 
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implicit in the 6j symbols.]] Instead, the conserved quantities are, as they must be, the total 
meson+baryon isospin and angular momentum, J tot and J to t- The S"-matrix element on 
the left-hand side is a physical partial-wave amplitude that can be compared directly with 
experiment. The "reduced S- matrix" under the summation is a meson- hedgehog amplitude, 
in slightly different notation than that of the previous Section. The relation between the 
two notations is: when the incoming and outgoing mesons are each pions, then Skkk — Sk, 
Sk,k-x,k-i = S^, Sk,k+i,k+i = S^, Sk,k-i,k+i = Sk,k+i,k-i = S{^>; when they are both 
cr's then Skkk = Sf 3 ; and when the incoming meson is a pion and the outgoing meson is a 
cr, then Sk,k-i,k = and Sk,k+i,k = S^; with all other elements vanishing. 



B. The Big- Small- Small- Big pattern 

For the remainder of this paper we specialize to the elastic case ttN — > ttN. For each value 
of L = L', there are then four a priori independent partial wave amplitudes, traditionally 
denoted L 2Itot 2J . For example, in the case of F-wave scattering (L = 3) the four physical 
amplitudes are F 15 , F 17 , F 35 , and F 37 . But to leading-order in large- N c , only two out of 
these four are independent. One can, for instance, solve for the two isospin-3/2 amplitudes 
as energy-independent linear combinations of the two isospin-1/2 amplitudes |^9|,^2j; this is 
an example of the It = Jt rule [ p6| , pT| . This holds in the Skyrme model, and because the 
group-theoretic expression (^) is the same, necessarily in the present <j-tt model as well. 
These relations are reasonably well obeyed by the experimental nN partial- wave data [22|, 
and are model-independent tests of large -/V c . 

Another interesting fact about the experimental data (see Fig. 4 ff., Ref. f22|): If for each 
L one juxtaposes the four amplitudes in the above order, namely L 12 l~i, Li,2L+i, -^3,2L-i 
and Lz^l+i, then they reveal a striking pattern termed the "Big-Small-Small-Big" pattern. 
Namely, the outer two amplitudes, Li j2 l-i and L 3>2 l+i, are characterized by relatively large 
excursions of the S-matrix element through the Argand circle, while the inner two am- 
plitudes, L12L+1 and L^ 2 l-i, show relatively much less motion. The Big-Small-Small-Big 
pattern is the single most consistent pattern characterizing the partial-wave S-matrix as a 
whole (the only clear exception to it being the D 35 ). 

Reproducing the Big-Small-Small-Big pattern is one of the noteworthy successes of the 



9 Note: if either the incoming or the outgoing meson is a a, then the associated 6j symbol has 
a zero in it and collapses to a product of Kronecker <5's. Conversely, the generalization of this 
expression to mesons that carry both isospin and spin, such as p's, involves 9j symbols, and is 
given in Ref. [24]. 
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Skyrme model | 2l| . It is equally well reproduced by the present <j-it model, as we illustrate 
in Fig. 17. In fact, the pattern emerges for the same dynamical reason |2~2 |: the fact that, 
for K > 1, in both the Skyrme model and in the a-ir model, the phase-shifts associated 
with Sfi" 1 are much smaller than those of S^ -1 (cf. Figs. 15-16). We therefore view it as a 
model-independent success of the large- N c approach, whether one chooses to use skyrmions 
or explicit baryon fields. 



C. The baryon spectrum of the large- N c a-rc model 

From the partial-wave amplitudes it is easy to extract the baryon resonance spectrum 
of the large-iV c <j-tt model. Rather than record when the phase-shifts cross 90 degrees (a 
crude criterion sensitive to background potentials), a more robust definition of a resonance, 
adopted by experimentalists, is to look for Lorentzian peaks in the "speed plots," i.e., the 
plots of \dTu tot j tot / du\ versus to. The speed plots for a few selected partial waves are depicted 
in Fig. 18. Some peaks are unambiguous, whereas others are admittedly "in the eye of the 
beholder," but the same can be said about the corresponding experimental data. 

Figure 19 displays the full resonance spectrum obtained in this fashion, through the 
if- waves (L = 5), limited to what we subjectively consider to be "two-star" resonances or 
better. The step-like structure, in blocks of alternating parity, is much more pronounced 
than in the Skyrme model, and certainly than in Nature. It can be partially accounted for 
by noting that, for L > 1, the reduced amplitudes of Fig. 16 dominate those of Figs. 14-15, 
so that for any fixed value of L, the resonance location in the four physical partial-wave 
amplitudes can be approximated by the resonance location in the single underlying reduced 
amplitude SfjT 1 . But this does not explain why the steps arrange themselves by definite 
parity (as we have indicated by the black bars below the horizontal axis), a feature for 
which we have no good understanding. 

In general, the resonances in the <j-tt model occur at substantially lower energies than in 
the Skyrme model, and in Nature. We have not explored the parameter space of our model 
[see Eq. (p8|)1 in an attempt to rectify this disparity (as we are confident could be done), not 
just because of the computationally- intensive character of these multi-stage calculations, but 
also due to the frankly "toy" intent of this model, which we have constructed for illustrative 
purposes. We are optimistic that a more realistic model, incorporating the vector mesons, 
and properly implementing chiral symmetry, would be in better agreement with the observed 
baryon spectrum, while posing no significant additional conceptual or numerical difficulties 
beyond those we have already confronted herein. 

The one parameter that we have experimented with is the nucleon size parameter a^, 
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defined in Eq. (|T^), which acts as an ultraviolet cutoff. A variation from our nominal value 
ax = 0.52 fm to = 0.60 fm shows no discernible effect on the resonance positions, and 
only slight changes in the Argand plots themselves, primarily in the P-waves, one of which 
is shown in Fig. 20. 



D. Some familiar problems 

We have seen that this large- N c a-n model (and, we presume, others like it with explicit 
baryon fields) shares some notable successes with the Skyrme model — the Big-Small-Small- 
Big pattern, the energy-independent relations between the J t ot = 1/2 and I tot = 3/2 nN 
amplitudes, the overall richness of the baryon resonance spectrum, etc. Moreover, the high- 
energy behavior of the partial wave amplitudes is much better than in the Skyrme model 
(see Sec. V). Not surprisingly, the a-ir model also shares some of the Skyrme model's failings. 
Figure 21 illustrates a specific partial wave amplitude in the a-7r model, juxtaposed with 
the experimental data. Obviously the real-world amplitude is much more inelastic than 
the present model. This is because, in the higher partial waves especially, multiple pion 
production soon dominates the experimental ttN amplitudes. Yet, formally, processes such 
as 7i N — > -kit-it N are down by powers of 1/N C compared with ttN — > ttN, and are therefore 
entirely absent from leading-order theoretical treatments such as the present paper — as 
well as from the leading-order skyrmion treatments [p)|-|2"3"|j , which share the same problem. 



Below the a threshold, the only source of inelasticity in the present model is the 7rA channel, 
exactly as in the Skyrme model. A theoretical means of summing at least some of the 1/N C 
corrections, namely those associated with multiple pion production, would immeasurably 
improve either approach. 

Just as serious is the failure of the o-tt model to bear even passing resemblance to 
experiment in the S and P waves. As is well known, these waves have been the 'Achilles 
heel" of the Skyrme model too. Interestingly, whereas the Skyrme model shows too few 
resonances in these waves, the o-ir model errs in the opposite direction: too many resonances, 
particularly in the P i3 and P 3i waves, and including spurious bound states in the and Su 



channels as already noted in Sec. V. The interested reader is referred to Ref. |22| for a lengthy 
discussion of the problems in these lower waves in the Skyrme model, which are related, in 
part, to the failure to incorporate the translational and (iso)rotational recoil of the hedgehog 
(formally 1/N C corrections, but numerically important). We expect that commentary to 
apply as well to models with explicit baryon fields. For example, the Weinberg- Tomozawa 



expression for the irN scattering lengths pEI] , which are predicted by current algebra, and 
which dominate the experimental S-w&ve amplitudes near threshold, formally appear only 
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at next-leading order in 1/N C |p2|j . This suggests that if one were to start from an improved 
effective hadron Lagrangian that respects chiral symmetry (we remind the reader that the 
present a-ir model does not), and if one were to calculate to next-leading order in 1/N C , the 
most glaring disagreement with experiment in the S'-waves ought to be repaired. Fixing the 
P-waves will require, at the least, (i) the splitting of the A from the nucleon (again, a 1/N C 
effect), and (ii) the incorporation of the Compton-type diagrams, particularly Figs, la and 
lb, the ameliorating effect of which has already been examined in the Skyrme model ]TT|,|T2" 
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APPENDIX A: NUMERICAL AND COMPUTATIONAL DETAILS 



1. Solving the coupled non-linear equations 

We describe here the method we use for solving the non-linear classical field equations, 
Eq. flilD, for the full a-ir model. (Exactly the same technique is used for the cr-only and tt- 
only warm-up problems discussed before that point.) These are coupled ordinary nonlinear 
differential equations. The nonlinearity implies that there is some sensitivity in finding 
solutions; indeed, for some ranges of parameters, one may not be able to find solutions at 
all. 

We are looking for solutions of these equations which are regular at the origin and which 
fall off asymptotically at large distances. Examining the indicial behavior of Eq. (f|l]) near 
r = 0, we find that 

F(0) = , F'(0) = A , G(0) = B G"(0) = . (Al) 

The indicial values A and B, along with asymptotic scale parameters C and D defined below 
in Eqs. (|A2|) and (|A"3|), are initially unknown. Their values will be fixed by the non-linearity 
when we solve the differential equations. 

As r — >• oo, the pion profile function is required to have the usual Yukawa-like fall-off, 

F(r) — ► Cexp(— m n r)/r . (A2) 

Things are more complicated for G(r) as r —>■ oo because the asymptotic behavior of G(r) 
is governed by the coupling term to two pions. That is, it will fall off like exp(— 2m n r)/r 2 
rather than the faster exp(— m a r)jr. After some analysis (solving a linearized Eq. Qf lb| ) 



using a Green's function technique with the F 2 {r) term on the right-hand-side providing the 
inhomogeneity) , 

G(r) -> D • - - £i[-(m CT + 2m>] , (A3) 

r Zrrirr r 



where Ei is an exponential integral p9| . 

The solution of Eq. ( ^T|) constitutes a two-boundary- value problem. We solve the equa- 
tions using a standard shoot and match technique PEfl . Not knowing, at first, the values 



for the "scale parameters" A. . . D, we make an initial guess for their values and proceed to 
refine them with an iterative procedure. 

The procedure is as follows. Given A . . . D, we integrate out from the origin using 
a Runge-Kutta technique to our matching radius, r m , which we choose to be ajv- We 
then Runge-Kutta integrate backwards to r m from r a , a point where the asymptotic forms 
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shown above are expected to hold. We typically choose r a to be 2 fm. The values at 
r m of F and G, and their derivatives, from the two integrations are then compared. The 
mismatches, or discontinuities, give four conditions which can be used to choose corrections 
to the guessed values of A ... D that would drive the discontinuities toward zero. (This is 
a generalization to four variables of the Newton-Rapheson method; it requires four more 
passes of Runge-Kutta integrations to compute the derivatives of the discontinuities with 
respect to the scale parameters.) Correcting the A . . . D as calculated, one can repeat this 
process, hopefully getting a better, less discontinuous F(r) and G{r). The procedure is 
iterated until it converges to a solution. 

The above iterative procedure was programmed in Fortran and was originally run as a 
batch job on a VAX- VMS minicomputer from a terminal command line. For small values of 
A (i.e., small non-linear contributions) the program would converge reasonably well. How- 
ever, as the non-linearity was increased, corresponding to values of A given by Eq. (p7|), 
the convergence became more delicate. Thus, we found it useful to make the code more 
interactive, so that the user could watch plots of F and G at every stage of the iteration 
and, if the process were going astray, stop it and start again with a new set of starting scale 
parameters. The computing was transferred to a NeXT workstation and a NeXTSTEP front 
end to the Fortran program was developed With this front end to aid the user, it was 
much easier finding a solution for A = 91. 



2. Runge-Kutta Approach to the Scattering Differential Equations 

The quantum scattering phase shifts in, say, the cr-only model of Sec. II, is given by the 
asymptotic form of the radial wave function determined by Eq. (p0|) . This linear differen- 
tial equation (and the coupled- channels variants of Sees. Ill and V) are readily solved by 
Runge-Kutta integration |j36| , simply by integrating out from the origin. Starting values 



are taken from the regular solution for the partial wave of the angular momentum K under 
consideration, i.e., with behavior like jxi.^) oc r K+1 . (The scale can be chosen arbitrarily 
because the equation is linear.) 

At a large enough distance, typically 10 fm or so, the solution is well-approximated by 
a linear combination of (<??") and n^(gr). Fitting the coefficients of these functions, one 
can then form the S- matrix element and thence compute the phase shift. For the coupled- 
channels case, the S'-matrix is found using the procedure given, for example, in Ref. |2T| 
utilizing matrices formed from these fitted coefficients. 
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3. Comments on Solving the Integral Scattering Equations 



Numerically, Eq. (|53|) is solved using the Nystrom method [36|. For our three-rowed 
column vector, this procedure involves the inversion of a large 3iV x 3iV matrix, where N 
is the number of meshpoints. The time to do that inversion goes like A 3 . We found that 
reasonable accuracy (1% or better) in the extracted S-matrix elements is obtained with 
N = 150 mesh points. These are typically distributed as follows: 120 mesh points, evenly 
spaced, from r = to 3 fm and 30 points, evenly spaced, from 3 to 10 fm. 

In debugging the coding we found it very useful to verify that the S'-matrix is not only 
unitary but also symmetric (a consequence of time reversal invariance). Another useful check 
on the code was to see that the scattering wave functions and S are continuous through the 
a threshold and agree with the prediction of the Runge-Kutta method above the threshold. 



APPENDIX B: FIGURE CAPTIONS 

1. Compton-type graphs contributing to meson-baryon scattering. Henceforth, directed 
lines stand for baryons, all other lines are mesons. The baryons B, B' and B" stand 
for three members of the I = J tower of baryons. Fig. 1(d) contains a purely mesonic 
loop and so is subleading in l/N c ; Figs. l(a)-(c) do not, and so contribute at leading 
order [9]. 

2. Exchange-type graphs contributing to meson-baryon scattering. Fig. 2(f) contains 
a purely mesonic loop and so is subleading in 1/N C ; Figs. 2(a)-(e) do not, and so 
contribute at leading order. In this paper we explicitly sum all such leading-order 
graphs in a model containing both tt and a mesons. 

3. Examples of exchange-type graphs contributing to (a) the baryon-baryon interaction, 
(b) the baryon-antibaryon interaction, and (c) the baryon-baryon-baryon interaction. 
As none of these particular graphs contains a purely mesonic loop, they all contribute 
to leading order in 1/N C , and can be summed using the methods of Ref. [9]. 

4. An uncrossed n-meson exchange graph included in the summation of Fig. 2. The 
shaded blobs contains meson self-interactions that do not concern us. 

5. A specific "crossing" or "tangling" of Fig. 4, for n — 4. There are n! such tanglings. 

6. The sum of all n\ tanglings including Figs. 4-5. In all three figures, the contents of the 
shaded blob are held constant, classical sources. 
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7. The tree-level one-point function a c \, where the baryon has been replaced by an exter- 
nal source j (x) . 

8. The fluctuating field <r qu propagating through the nontrivial background generated by 
(T c i. Vertices can be read off from the quadraticized Lagrangian, Eq. (15). 

9. The profile function <7 c i(r) = G(r) for the cr-only model of Sec. II. 

10. Phase shifts for the <7-only model, as a function of a momentum, for model parameters 
m a = 600 MeV, g = 13.6, and = 0.5 fm. 

11. Classical profile F(r) in inverse fermis for the 7r-only model of Sec. III. 

12. Classical profile functions F(r) and G(r) for the a-n model, with parameters as in 
text. 

13. Diagonal and off-diagonal potentials appearing in the small-fluctuations equations. 

14. The phase shifts 5k plotted against meson momentum k in the baryon rest frame 
(recall that baryon recoil is subleading in 1/N C ). 

15. The phase shifts 5^ plotted against meson momentum k in the baryon rest frame. 

16. The phase shifts 5% 2 plotted against meson momentum k in the baryon rest frame. 

17. Argand plots for the four F-wave elastic ttN partial wave amplitudes, illustrating the 
"Big-Small-Small-Big" pattern (the Fi 5 and F 37 are big, the others small). 

18. Examples of speed plots for selected irN amplitudes, from which resonance positions 
are extracted. 

19. The baryon spectrum of the a-ir model. The vertical axis measures excitation energy 
in MeV above the nucleon/hedgehog mass. Not pictured are the two spurious bound 
states in the Su and S31 channels. 

20. Dependence of the Pu amplitude on the nucleon size parameter ajy, for 0.52 and 0.60 
fm (solid and dashed lines, respectively). 

21. Comparison of Gn amplitude between the a-n model (solid line) and experiment 
(dotted line). The pion kinetic energy for each curve ranges from to 1600 MeV. 
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FIG. 1. Compton-type graphs contributing to meson-baryon scattering. Henceforth, directed 
lines stand for baryons, all other lines are mesons. The baryons B, B' and B" stand for three 
members of the I = J tower of baryons. Fig. 1(d) contains a purely mesonic loop and so is 
subleading in 1/JV C ; Figs. l(a)-(c) do not, and so contribute at leading order [9]. 
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FIG. 2. Exchange-type graphs contributing to meson-baryon scattering. Fig. 2(f) contains a 
purely mesonic loop and so is subleading in 1/N C ; Figs. 2(a)-(e) do not, and so contribute at leading 
order. In this paper we explicitly sum all such leading-order graphs in a model containing both tt 
and o mesons. 
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FIG. 3. Examples of exchange-type graphs contributing to (a) the baryon-baryon interaction, 
(b) the baryon-antibaryon interaction, and (c) the baryon-baryon-baryon interaction. As none of 
these particular graphs contains a purely mesonic loop, they all contribute to leading order in 1 /N c , 
and can be summed using the methods of Ref. [9]. 




FIG. 4. An uncrossed n-meson exchange graph included in the summation of Fig. 2. The 
shaded blobs contains meson self-interactions that do not concern us. 
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FIG. 5. A specific "crossing" or "tangling" of Fig. 4, for n = 4. There are n! such tanglings. 
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FIG. 6. The sum of all n\ tanglings including Figs. 4-5. In all three figures, the contents of the 
shaded blob are held constant, classical sources. 
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FIG. 7. The treedevel one-point function a c \, where the baryon has been replaced by an external 
source j(x). 
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FIG. 8. The fluctuating field a„ n propagating through the nontrivial background generated by 
a c \. Vertices can be read off from the quadraticized Lagrangian, Eq. (15). 
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FIG. 9. The profile function <J c \{r) = G(r) for the cr-only model of Sec. II. 
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FIG. 10. Phase shifts for the cr-only model, as a function of a momentum, for model parameters 
m a = 600 MeV, g = 13.6, and ajy — 0.5 fm. 
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FIG. 11. Classical profile F(r) in inverse fermis for the 7r-only model of Sec. III. 
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FIG. 12. Classical profile functions F(r) and G(r) for the a-it model, with parameters as in 
text. 




FIG. 13. Diagonal and off-diagonal potentials appearing in the small-fluctuations equations. 
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FIG. 14. The phase shifts Sx plotted against meson momentum k in the baryon rest frame 
(recall that baryon recoil is subleading in 1/N C ). 





FIG. 17. Argand plots for the four F-wave elastic ttN partial wave amplitudes, illustrating th 
Big-Small-Small-Big' 1 pattern (the F15 and F37 are big, the others small). 
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FIG. 18. Examples of speed plots for selected irN amplitudes, from which resonance positions 



are extracted. 
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FIG. 19. The baryon spectrum of the a-ir model. The vertical axis measures excitation energy 
in MeV above the nucleon/hedgehog mass. Not pictured are the two spurious bound states in the 
Sn and S31 channels. 
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FIG. 20. Dependence of the P\\ amplitude on the nucleon size parameter u \. for 0.52 and 0.60 
fm (solid and dashed lines, respectively). 
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FIG. 21. Comparison of Gn amplitude between the <j-tt model (solid line) and experiment 
(dotted line). The pion kinetic energy for each curve ranges from to 1600 MeV. 
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